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We present the quantum phase transition in two capacitively coupled arrays of superconducting 
quantum dots (SQD). We consider the presence of gate voltage in each superconducting island. 
We show explicitly that the co-tunneling process involves with two coupled SQD arrays, near 
the maximum charge frustration line is not sufficient to explain the correct quantum phases with 
physically consistent phase boundaries. We consider another extra co-tunneling process along each 
chain to explain the correct quantum phases with physically consistent phase boundaries. There 
is no evidence of supersolid phase in our study. We use Bethe-ansatz and Abelian bosonization 
method to solve the problem. 
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Josephson junction arrays have attracted considerable 
interest in the recent years owing to their interesting 
physical properties [l-14|. Currently such arrays can 
be fabricated in restricted geometries both in one and 
two dimensions with well controlled parameters 
At the same time this is one of the paradigms to study 
the physics of quantum phase transition Hil. The 
Coulomb charging energy of the system is due to the 
small capacitance of the grains or junction that domi- 
nate strong quantum phase fluctuations of the supercon- 
ducting order parameter and may drive the system into 
the insulating state at zero or very low temperature. The 
quantum fluctuations are controlled by the parameters of 
the system such as charging energy of the superconduct- 
ing quantum dots (SQD) and Josephson coupling {Ej ) 
between them. 

In two dimensions, the universality class of these tran- 
sitions has already been investigated in detail lol- 14|. 
Nevertheless there are also quantum phase transition sce- 
nario that can take place in a Josephson junction ladder 
system as shown in the Fig. 1, where we consider two ca- 
pacitively coupled one dimensional SQD arrays. In this 
simplest two-leg ladder system, we will be able to pre- 
dict different quantum phases due to the interplay be- 
tween Coulomb charging energy and Ej. Here we also 
show explicitly the role of co-tunneling effect. For large 
Coulomb charging energy sequential tunneling is not en- 
ergetically favored process in the system. The major tun- 
neling process in the system occurs via the co-tunneling 
process [l9| . Its origin is quantum mechanical. In this 
process, tunneling occurs through the virtual state with 
energy equal to the on-site Coulomb charging energy. 
An appropriate co-tunneling process extract the correct 
quantum phases in the system. Otherwise, it leads to 
the wrong analysis of the system as we have found in 
Ref . [20I . In our previous studies we only emphasis the 
one dimensional SQD array, we only explain the impor- 
tance of next-nearest-neig hbor (NNN) i;j 0, II . 



In this paper we study the quantum phase transition of 
capacitively coupled two arrays of SQD and they are con- 
nected through Josephson junctions. We consider the 
presence of gate voltage in each SQD which introduce 
the charge frustration in the system. At the charge frus- 
tration line the Coulomb charging energy of the system 
is degenerate for the difference of one Cooper pair in the 
island. Its also refer in the literature as a charge degener- 
acy point. Our main motivation is to study the different 
quantum phases around the maximum charge frustration 
line. Choi et al. [2^ have also studied the same problem 
They have shown that in a coupled chain the major trans- 
port along both the chains occurs via co-tunneling of the 
electron-hole pairs They have not done any rigorous ana- 
lytical exercise based on techniques, which are applicable 
in low dimensional quantum many particle systems and 
that seed mistakes. There are three major mistakes in 
that paper (1) There are two kinds of Luttinger liquid 
(LL) phase (describe as RLl and RL2) but they have ob- 
tained only one. (2) Incorrect quantum phase analysis 
and physically inconsistent phase boundaries due to lack 
of correct analytical derivations and physical interpreta- 
tions. (3). There is no evidence of super solid (SS) phase 
in this model. 

We show explicitly in this study that a single co-tunneling 
process in the two leg SQD ladder system is not sufficient 
to produce the correct quantum phases with physically 
consistent phase boundaries. We also show that the co- 
tunneling process along each chain is also necessary for 
the correct phase boundaries. In our model system, su- 
perconducting islands are connected with Ej. The charg- 
ing energies along each array are the on-site Coulomb 
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charging energy due to the self capacitance, Eq — 

2 

and the junction charging energy, Ei = ^r- due to junc- 
tion capacitance Ci . The arrays are coupled to each other 
via the capacitance C/, it contributes charging energy, 
Ej — ^Tj- . But there are no Josephson coupling between 
the two arrays. The array capacitances are small, which 
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FIG. 1: This is the schematic diagram of our model systems. 
Two arrays of superconducting quantum dots are connected 
through the capacitance, Ci. Co is the self capacitance of the 
dot and Ci is the junction capacitance. In our superconduct- 
ing circuit, small circles denote the dots and crosses denote 



the Josephson junctions, 
each dot. 



Vg is the gate voltage applied in 



leads to Ej << Eo,Ei. Therefore each arrays are in 
the insulating phase in the absence of co-tunneling effect. 
Here we consider the limit Ej « Eo,Ei. The applied 
gate voltage induce a charge in every superconducting is- 
land by an amount, Ug = '^^J" ; this applied gate voltage 
introduce the charge frustration in the system. We do 
the quantum phase analysis near the maximum charge 
frustration line (1%— TV — 1/2| << 1), around this line 
we obtain several interesting quantum phases which have 
not been noticed in the previous studies d, . The 

Hamiltonian of the system is given by 

H = 2e^ ^ {ni{x) - Ug) Cj^^]{x,x ) {ni{x ) - Ug) 

ll' ^XX 

-EjY,cos{di{x)~Qi(x^\)) (1) 

ni (x) is the number of Cooper pairs and 0/ (x) is the phase 
of the superconducting order parameter at the site x of 
the 1th array. ni(x) and di{x) are the quantum mechani- 
cally conjugate variables. One can write the capacitance 
matrix in the block form. 

C~l, ix,x' )^ Cix, X )A + 5^^^, CiB (2) 



Where 



A 



B = 



1 
1 



1 -1 
-1 1 



C{x,x) = CoS^^^. + Ci[25,y-<5,y+i-^_,_,] (3) 



(4) 



with the components 

i?c<°^ =(7o^[n+(:r)-2n,]VFo5^(n-(a;)f (5) 

X X 

Hc^'^ = UiY.[n+{x)-2ng][n+{x + l)-2ng] 

X 

+Vi^n^{x)n^{x + 1) (6) 



Hj = ~EjJ2cos{0i{x) -ei{x + 1)) 



(7) 



Where n±{x) = ni{x) ± ^2(2;) and the coupling 



strengths are given hy U ^ 2Eq, Ui 



'Co 



E, 



and Vi - ^=^. The first term of the Hamiltonian iJc'"' 
contains an important message regarding the difference 
between the charging energy at near to the maximum 
charge frustration line and the particle-hole symmetric 
line {ug = 0). The charge configurations which do not 
satisfy the condition for n+(x) = 1 at the maximum 
charge frustration line generates a gap in the excitation 
spectrum of the order of on site charging energy. The 
ground state of iJc^°'' separated from the excited state 
by the gap of the order of Ei. This excited state has two 
fold degeneracy for every values of x, corresponding to 
= ±1. It is convenient to work in the charge con- 
figuration with = and n_(a;) = ±1 is termed as 
a reduce Hilbert space of the problem. Presence of finite 
Ej lifts this degeneracy and the ground state of Hc^ is 
mixed with the state n_(a:) = ±2. Now the relevant re- 
duce Hilbert space, n+(a;) — and n^{x) — 0, ±2. Here 
we mention the effective model of our system, wherein 
we follow the references . They have found the effective 
Hamiltonian up to the second order in 



eff 



= P[H + Hj 



1 -P 



oHj]P 



(8) 



E - He" 

Here P is the projection operator onto the reduce space 
and finally they have obtained the effective Hamiltonian. 



Heff = -,J}_^S^-{x)S%x + l) 

X 

-J 



(9) 



where the XY component of exchange interaction and 
the z component of exchange anisotropy are respectively 
J — and 7 — ■ No analytical expression 

for is given in Ref. . Here we estimate — §rMr-. 

. . . Cl-tyJ 

Pseudo spin operators of the above Hamiltonian are the 
following: 



s^{x) = p!^M-I^P 

S+{x) = Pe-'^i(^)(l -F)e-'^^(^)P 
S-{x) = Pe-'^^(^)(1 -P)e-^^i(^)p. 



(10) 
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Now we would like to explain the pseudo-spin terms 
of the Hamiltonian in terms of charge representation 
from where we start. The first term of the Hamiltonian 
originates from the first term of Eq. 9 which the 
Coulomb charging energy of the SQD in the form of 
n{x). The second term of the Hamiltonian originates 
from the co-tunneling process, i.e, from the second term 
of Eq. 9. They have obtained the Hamiltonian in an 
effective one dimensional Hamiltonian but they have not 
done any correct analytical calculations based on the 
Abelian bosonization and Bethe ansatz methods which 
are suitable for one dimensional quantum many body 
system. Here we would like to do the quantum phase 
boundaries analysis based on the exact Bethe ansatz 
calculation based results. We will present the short 
comings of this effective Hamiltonian from that analysis. 



S^{x) = - 1/2 , 

(x) — "ipix) exp[i7r rij] , 

j=-oo 

S+{x) ^ ip\x) exp[-iTT uj] , (11) 



K=l/4 



K=l 



j = -oo 



where n{x) = {x)ip{x) is the fermion number at the 
site x. 



(7/;t(a;-Hl)V(x-M)-l/2) 

^ Y (V-^x + l)^(x) -I- h.c.) (12) 



In order to study the continuum field theory of these 
Hamiltonians, we recast the spinless fermion opera- 



tors in terms of field operators by this relation [21 1 



■>p{x) 



ikpx ^^(x)] where ipB.{x) 



and ipLix) describe the second-quantized fields of right- 
and the left-moving fermions respectively. We want to 
express the fermionic fields in terms of bosonic fields by 
the relation ibJx) = e"* ('■'^("^^ - ''(^^ , where r 

is denoting the chirality of the fermionic fields, right (1) 
or left movers (-1). The operators commute with the 
bosonic field. Uj. of different species commute and Ur of 
the same species anti commute. (f> field corresponds to the 
quantum (bosonic) fiuctuations of spin and 9 is the dual 
field of 0. They are related by the relation (pR, ~ 6 — cj) 
and (pL — + (j). After doing the continuum field 
theory exercise, Heff become 

H,ff ^Ho + ^ J dx : cosiA^i(bix)) : (13) 




K=l/4 K=l 



FIG. 2: We present the one dimensional schematic diagram 
for quantum phases as a function of Luttinger liquid param- 
eters (K). In the figure, we present also values of K at the 
phase boundaries. The value of K at the wings of the arrows 
denote the value of K at the phase boundaries. CDW is the 
charge density wave phase. RL2 is the repulsive Luttinger liq- 
uid phase of second kind. SC is the superconducting phase. 
The limit K — 2 only achieve when the system is in the first 
order commensurability, i.e., each dots are occupied with inte- 
ger number of Cooper pairs. MI is the Mott-insulating phase 
of the system. The presence of gate voltage is not explicit in 
K, therefore the RLl phase not appears in figure, we can only 
achieve this phase when the applied gate voltage exceed the 
CDW gap. 



Ho = ^ J dx[{dxd)'^ + {dx<p)'^] is the non-interacting part 
of the Hamiltonian. The collective velocity of the system 
(v) and Ki are the two LL parameter. We use exact 
Bethe ansatz solution to calculate 



TT 4- 2sin -'-7 



(14) 



For J < with relatively small applied gate volt- 
age, the antiferromagnetic Ising interaction dominates 
the physics of anisotropic Heisenberg chain. When the 
applied gate voltage is large the chain is in the ferromag- 
netic state. In the language of interacting Cooper pairs, 
the Neel phase is the Charge density wave phase with 
period 2 i.e there is only one Cooper pair in every two 
sites. For the ferromagnetic phase system is in the Mott 
insulating phase. The emergence of two LL phases can be 
ascribed to the following reasons: RLl and RL2, respec- 
tively occur due to the commensurate-incommensurate 
transition and the criticality of the Heisenberg XY model. 
The system is either in RLl for if < 1/2 or in RL2 for 
K > 1/2. The physical significance of RLl is that the 
coupling term is relevant but the applied gate voltage 
on SC island breaks the CDW gapped phase, whereas in 
RL2, system is gapless for Ki > 1/2. 
In Fig. 2, we present quantum phases of our study as a 
function of K. We apply the Luther-Emery (LE) trick 
22] in the massive phase of sine- Gordon field theory to 
evaluate Ki (=1/4) at the phase boundary between the 
CDW and RLl and also for CDW and RL2. RLl phase 
is not explicit in the figure because we only achieve this 
phase under the application of gate voltage when it ex- 
ceed the gap of CDW state. The value of Ki at the 
phase boundary between RL2 and CDW phase is 1/4. 
From the analysis of Ki we obtain IQX^Ei^ = —Ej^, this 
condition is unphysical because A, Ej and Ej are (+ve) 
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quantities. The value of isT is 1 at the phase boundary be- 
tween the RL2 and superconducting phase. The analysis 
yields the condition 16A^i?/^ = 0, which is again unphys- 
ical. There is no chance to get superconductivity in this 
model Hamiltonian. We surprise that without doing any 
correct quantum analysis of this model Hamiltonian, the 
authors of Ref. (2^ have claimed the existence of SC. 
To get the correct physical behavior between the differ- 
ent quantum phases, one has to be considered the Co- 
tunneling processes along each chain separately. We have 
realized during our calculations that to get an attrac- 
tive interaction between the Jordan- Wigner (spinless) 
fermions, we have to consider the higher order expan- 
sion in This higher order expansion leads to the 
virtual state with energies exceeding Eco- In this second 
order process, extra contribution appears as 



Hp 



-3E, 



4Ei 

X 

-^J2(^H^ + 2)S-{x) + h.c) (15) 
J^o — 

X 

[23j . The total effective Hamiltonian of the system 
under the combined co-tunneling process is 

-^J2(^Hx)S^{x + i) + h.c) 

X 

-^Y.^S\x + 2)S-{x)+h.c) (16) 

^0 X 

After doing the quantum field theory, we get the effective 
Hamiltonian 



Heff - Hq- 



7 J 



3Ej 



4Eo 



27ra2 



dx : cos{4^/k^(I){x)) : (17) 



LL parameter of Hef / is 



K2 = 



J 



J + 4:/tt{jJ 



3E., 



(18) 



4Eo 



As we notice from the Fig. 2 that value of K2 = 1/4 
at the phase boundary between RL2 and CDW phase. 
The parametric condition at the phase boundary is the 
following ~ + 3 This condition is consistent 

physically. Therefore we prove the importance of this 
extra co-tunneling process. We see from the Fig. 2 that 
value of the is 1 at the phase boundary between the 
RL2 and superconducting phase. The parametric anal- 
ysis for this value of K2 yields \%}?Ei^ — "iEj^ , we 
obtain, Ej^ = %\^Ei'^. This is a physically realizable 
condition for phase boundaries. It is clear from the Fig. 
2 that there is no simultaneous presence of CDW phase 



and SC phase. SC phase occurs only when K > \ and 
CDW phase occurs only when 1/2 < if < 1. Therefore 
there is no evidence of SS phase for this model system. 
Conclusions: We have found the all correct quantum 
phases of this model Hamiltonian based on the rigorous 
analytical calculations. We obtain the condition for phys- 
ically consistent phase boundaries which were absent in 
the previous studies. There is no evidence of super solid 
phase. 
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We follow the Luther-Emery trick during 
the analysis. One can write the sine-Gordon 
Hamiltonian for arbitrary commensurability as 
H2 = Hoi + X J dx cos{2n \fK 0(2:)), where 
n is the commensurability and A is the coupling 
strength, i/01 is the free part of the Hamilto- 
nian. We know that for the spinless fermions, 
tpR^'i'L + ipL^'i'R = / dxcos(2y/K(l){x)), which 

is similar to the analytical expression of sine-Gordon 
coupling term but with the wrong coefficient inside the 
cosine. One can set 4>{^) = 2\/l{(j){x) then the above 
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equation become i?3 = Hoi + X J dxcos{24>{x)). K and [23] D. Giuliano, and P. Sodano, Nucl. Phys. B 711, 480 
K are related by the relation, K = ^. At the phase (2005). 
boundary, K = 1 that implies K = l/n^. 



